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Abstract 
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■ The nuclear density functional framework, based on chiral dynamics and the symme- 

. try breaking pattern of low-energy QCD, is extended to the description of collective 

nuclear excitations. Starting from the relativistic point-coupling Lagrangian intro- 
duced in Ref. [1], the proton-neutron (quasiparticle) random phase approximation 
is formulated and applied to investigate the role of chiral pion-nucleon dynamics in 
excitation modes involving spin and isospin degrees of freedom, e.g. isobaric analog 
' states and Gamow- Teller resonances. 

O ' PACS: 21.30.Fe, 21.60.Jz, 24.30.Cz 
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3 1 1 Introduction. 

> 

^ ■ In a series of recent articles [1,2,3,4,5,6] we have developed a novel approach 

to nuclear density functional theory which emphasizes the relationship be- 
tween aspects of low-energy, non-perturbative QCD and the nuclear many- 
body dynamics. A relativistic nuclear energy density functional has been for- 
mulated, starting from two basic features which establish the link with low- 
energy QCD and its symmetry-breaking pattern: (a) strong scalar and vector 
fields related to in-medium changes of QCD vacuum condensates; (b) long- 
and intermediate-range interactions generated by one-and two-pion exchange, 
derived from in-medium chiral perturbation theory, with explicit inclusion of 
A (1232) excitations and in combination with three-nucleon (3N) interactions. 

The resulting nuclear matter energy density functional [6] is mapped [3] onto 
the exchange-correlation energy density functional of a covariant point-coupling 
model for finite nuclei, including second order gradient corrections. In this 
approach the nuclear ground state is determined by the self-consistent solu- 
tion of the relativistic generalizations of the linear single-nucleon Kohn-Sham 
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equations. The construction of the density functional involves an expansion of 
nucleon self-energies in powers of the Fermi momentum up to and including 
terms of order k®. In comparison with purely phenomeno logical mean-field 
models, the built-in QCD constraints and the explicit treatment of pion ex- 
change restrict the freedom in adjusting parameters and functional forms of 
the density-dependent couplings. This approach to nuclear dynamics has been 
tested in an analysis of ground-state properties of a broad range of spherical 
and deformed nuclei [1]. The results have been compared with experimental 
data on binding energies, charge radii, neutron radii and deformation param- 
eters for several isotopic chains. It has been demonstrated that the results for 
nuclear ground-state properties are at the same level of quantitative compar- 
ison with data as the best phenomenological relativistic mean-field models. 

In this work the low-energy QCD-constrained nuclear density functional frame- 
work will be extended to the description of excited states and, in particular, to 
spin-isospin collective excitations. Starting from the relativistic point-coupling 
Hartree-Bogoliubov model of Ref. [1], we formulate the proton- neutron (quasi- 
particle) random phase approximation (PN-QRPA). This is then applied in 
an analysis of isobaric analog states and Gamow- Teller resonances in doubly- 
closed and open-shell nuclei. Our aim is to study the role of chiral pion-nucleon 
dynamics in charge-exchange modes. The results will be compared with avail- 
able data, and with QRPA calculations performed using one of the modern 
phenomenological relativistic density-dependent effective interactions based 
on the meson-exchange representation of effective nuclear forces. 



2 The Quasiparticle Random Phase Approximation 

The matrix equations of the relativistic random phase approximation (RRPA) 
for an effective Lagrangian with explicit density dependence of the meson- 
nucleon couplings have been derived in Ref. [7], and for nucleon-nucleon point- 
coupling effective Lagrangians in Ref. [8]. The relativistic QRPA for collec- 
tive excitations in open-shell nuclei has been formulated in Ref. [9] using the 
canonical single-nucleon basis of the relativistic Hartree-Bogoliubov (RHB) 
model. The extension to spin and isospin excitations (PN-QRPA) has been 
described in Ref. [10]. Our implementation of the PN-QRPA is identical to 
that of Ref. [10], with the important difference that it is now based on the 
point-coupling effective Lagrangian with density-dependent couplings deter- 
mined by chiral pion-nucleon dynamics [1], rather than on phenomenological 
meson-exchange effective interactions. Here we only collect the essential ex- 
pressions and refer the reader to Refs. [7,8,9,10] for technical details of the 
relativistic QRPA. 

We consider transitions between the + ground state of a spherical even-even 
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parent nucleus and a state with angular momentum and parity J 77 of the odd- 
odd daughter nucleus. The spin-isospin dependent coupling terms are gener- 
ated by the following interaction parts of the Lagrangian density r~l : 



£int — £-IV + £nN + £-LM (1) 

£iv = -\(%l{p)<&N>) ■ <?m) - \p { ^{p)i^h^) ■ (2) 
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£nN = -—^1^ ■ (3) 



C-lm = 2^o(p)(^^757m^) • (^75 W) • (4) 

The first term Crv collects the isovector parts of the nuclear density func- 
tional of Ref . [1] , with the density-dependent coupling strengths G^s an d G^ v 
completely determined by chiral pion-nucleon dynamics. Details are given in 
Ref. [1]. The same approach, without further adjustment, will also be used for 
the proton-neutron particle-hole (ph) residual interaction. Because of parity 
conservation, the one-pion direct contribution vanishes at the mean-field level 
in the calculation of a nuclear ground state. The pion-nucleon Lagrangian 
density C n N Eq. (3) must be included, however, in calculations of excitations 
that involve spin and isospin degrees of freedom [1 1] . The axial vector coupling 
constant: qa = 1.3, and the pion decay constant: f n = 92.4 MeV. 



Short-distance spin-isospin dynamics is encoded in the zero range Landau- 
Migdal term Clm, Eq. (4). In the non-relativistic limit the corresponding 
two-body interaction reduces to the familiar form g' cri ■ (t 2 ti - f^. In a recent 
detailed study of the role of 27r-exchange in the interaction of quasi-nucleons 
at the Fermi surface \pi\ = \f>2\ = kf [12], the four Landau-Migdal parameters 
fo(kf), fd(kf), go(kf) and g' (kf) which characterize the four spin-isospin chan- 
nels of the isotropic (I = 0) part of the N-N quasiparticle interaction, have 
been calculated using in-medium chiral perturbation theory (ChPT). This cal- 
culation includes contributions from l7r-exchange, iterated l7r-exchange and 
27r-exchange with virtual A(1232)-isobar excitations. The resulting depen- 
dence on the Fermi momentum kf is converted into density dependence using 
the relation for symmetric nuclear matter: p = 2fc^/(37r 2 ). 

In general, the explicit density dependence of the couplings introduces addi- 
tional rearrangement terms in the residual two-body interaction of the RRPA 
[7]. However, since these terms include the corresponding isoscalar ground- 
state densities, it is easy to see that they are absent in the charge-exchange 
channels, and the residual two-body interaction reads: 



Vectors in isospin space are denoted by arrows, and boldface symbols indicate 
vectors in ordinary three-dimensional space. 
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V ph ( ri , r 2 ) = GP s {p) fy ■ f 2 (l)i(l)2 5 3 (n - r 2 ) 

+ 4t(p) n • r 2 (/3 7 M )i(/^) 2 5 3 (n " r 2 ) 

n 2 p -m w \r 1 -r2\ 

+ Tfif2 * ■ (Si • V)(S 2 • V)- 

16tt/2 \r 1 -r 2 \ 

+ 0o(p) n ■ r 2 (^E)i ■ (/?£) 2 5 3 ( n - r 2 ) . (5) 

Here S is the 4x4 matrix representation of the spin operator. We note 
that there is no double counting of the one-pion contribution, because the 
Hartree diagram does not contribute to the isotropic part of the quasi-nucleon 
interaction in the long wavelength limit, i.e. to g' Q (p) [12]. 

Pairing interactions in the particle-particle channel have also been derived 
using in-medium chiral dynamics [13]. However, for the convenience of quanti- 
tative comparison with previous relativistic PN-QRPA calculations for open- 
shell nuclei of Ref. [10], and consistent with the model that we have developed 
for ground-state properties in Refs. [1,2,3], in the T = 1 channel the pairing 
part of the Gogny force is used 

VZlin, r 2 ) = £ e -«n-r2)/^f ( W . + B .pc _ H .pr _ Mi P°P T ), (6) 
i=l,2 



with the set D1S [14] for the parameters /ij, Wi, B i: Hi and Mj (i = 1,2), 
and a short-range repulsive Gaussian combined with a weaker longer-range 
attractive Gaussian in the T = channel: 

V^o(ri, r 2 ) = -V £ 9j e~^-^ 2 H 5=1 , T=0 , (7) 



where tl s=1T=0 projects onto states with S = 1 and T = 0. The ranges of the 
two Gaussians /ii=1.2fm and /i 2 =0.7fm are the same as for the Gogny inter- 
action Eq. (6), and the choice of the relative strengths g\ = 1 and g 2 = —2 
makes the force repulsive at small distances. The overall strength parameter 
of the T = pairing interaction Vq = 250 MeV is the same as in the rela- 
tivistic PN-QRPA calculation of Gamow- Teller resonances of Ref. [10]. The 
T = 1 Gogny interaction Eq. (6) is also used in the pairing channel of the 
RHB equations which determine the ground state of the initial nucleus. The 
phenomenological particle-particle interaction used here differs only slightly 
from the one derived from in-medium ChPT [13]. A complete investigation of 
in-medium pionic fluctuations (in the particle-hole and in the particle-particle 
channel) will be presented in a forthcoming paper. 

The two-quasiparticle configuration space includes states with both nucleons 
in discrete bound levels, states with one nucleon in a bound level and the 
second nucleon in the continuum, and also states with both nucleons in the 
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continuum. In addition to configurations built from two-quasiparticle states of 
positive energy, the RQRPA configuration space includes pair-configurations 
formed from the fully or partially occupied states of positive energy and the 
empty negative-energy states from the Dirac sea. In Refs. [9,10] it has been 
shown that the inclusion of configurations built from occupied positive-energy 
states and empty negative-energy states is essential for the consistency of 
the relativistic (proton-neutron) QRPA (current conservation, decoupling of 
spurious states, sum rules). 

The total strength for the transition between the ground state of an even-even 
spherical (N,Z) nucleus and the excited state \XJ) of the corresponding odd- 
odd (N+1,Z-1) or (N-1,Z+1) nucleus, induced by an isovector single-particle 
operator T JM , reads 



E 

pn 



(p\\T J \\n) 



7 Y XJ v u 

1 pn v p u n 



where p and n denote proton and neutron quasiparticle canonical states, 
respectively, with occupation factors u p , v p , u n , v n . X XJ and Y XJ are the 
forward- and backward-going QRPA amplitudes for the state |AJ). The dis- 
crete strength distribution is folded with the Lorentzian function 



In the calculations performed in this work the choice for the width of the 
Lorentzian is 1 MeV. 



2.1 Determination of g' (kf). 



The density-dependent Landau-Migdal parameter g' (kf) of Ref. [12] can be 
expanded in powers of the Fermi momentum kf to the same order as G^g(kf) 
and G$r(k f ) [1]: 

^(M = ci + c 2 ^ + c 3 ^ 2 + c 4 ^y , (10) 

with the coefficients q listed in Tab. 1, and A = 27rf n ~ 0.58 GeV chosen as 
a reference scale. We note that, in addition to the terms which are completely 
determined by the in-medium ChPT calculation of pion-exchange diagrams, 
the constants c%, C3 and C4 involve three parameters as displayed in Tab. 1: 63, 
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Table 1 

Coefficients in the expansion of the spin-isospin Landau-Migdal parameter g'o(kf) 
Eq. (10). The best-fit constants used in the present calculation (middle column), 
are shown in comparison with the estimates of Refs. [6,12] (right column). 



coefficients 


constants 


Cl 


(fm 2 ) 


-1.80 fm 2 - £b 3 




-2.93 




3.05 [6] 


C2 


(fm 2 ) 


1.17 fm 2 










C3 


(fm 2 ) 


-9.84 fm 2 + j^bl T 


w = 


-7.30 




8.03 [12] 


C 4 


(fm 2 ) 


21.13 fm 2 - £b% r 




-8.45 




B.45 [12] 



65 T and feg r , respectively, which encode the additional short-range dynamics 
of the NN and 3N interactions. In Ref. [6] the subtraction constant 63 = 
—3.05 was adjusted to empirical nuclear matter properties at saturation point. 
Taking into account ground-state properties of finite nuclei in addition, this 
constant has been fine-tuned to 63 = —2.93 in Ref. [1]. This latter value will 
consistently be used both in the Dirac Hamiltonian which determines the 
nuclear ground state, and in the isovector residual RPA interaction Eq. (5). 
The remaining two short-distance constants b^ T and 6g T cannot be constrained 
by ground-state properties of spin-saturated nuclear matter. In Ref. [12] the 
value of b^ T has been adjusted so that the resulting g' is consistent with 
empirical values at saturation density. Here we retain the value 6g r = —8.45 
from Ref. [12], and adjust b^ T in such a way that the PN-RPA calculation 
reproduces the excitation energy of the Gamow- Teller giant resonance (GTR) 
in Pb. In Fig. 1 we plot the Landau-Migdal parameter function of 

the Fermi momentum kf for several values of b^ T in the interval [—5.1, —9.5]. 
In the panel on the right the corresponding Gamow- Teller strength functions 
in 208 Pb. The experimental position of the GTR in 208 Pb: 19.2 MeV [15,16,17], 
is reproduced with b^ T = —7.30. This value is remarkably close to the estimate 
of Ref. [12]. 

With the constants b 3 , b^ T and 6g r adjusted as described above, we plot the 
resulting; function of the Fermi momentum in Fig. 2. The values around 

saturation point (kf ~ 260 MeV) are compared with: (a) a Brueckner calcula- 
tion with a phenomenological iViV-interaction [18] (dark gray rectangle); (b) 
a recent estimate based on the universal low-momentum potential Vi ow ^k [19] 
(black rectangle); (c) an empirical estimate based on the quenching of the 
Gamow- Teller strength [20] (light gray rectangle). We notice good agreement 
between the predicted values of g' Q at saturation density, even though they 
have been obtained with very different models or estimated empirically. The 
functional dependence for kf > 200 MeV is in qualitative agreement with the 
Brueckner-Hartree-Fock calculation of Ref. [21] (black dots), which has em- 
ployed the Argonne Vi 8 two-body potential [22], and a three-body force from 
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Ref. [23]. 



3 Applications: spin-isospin collective excitations 

Collective spin and isospin excitations in nuclei have been the subject of nu- 
merous experimental and theoretical studies (for extensive reviews see Refs. [24,25]). 
Nucleons with spin up and spin down can oscillate either in phase (spin scalar 
S=0 mode) or out of phase (spin vector S=l mode). The spin vector, or spin- 
flip excitations can be of isoscalar (S=l, T=0) or isovector (S=l, T=l) nature. 
These modes provide information on the spin and spin-isospin channels of ef- 
fective nuclear interactions. 

3.1 Isobaric analog resonances 

As a first application of this new PN-QRPA model, we calculate the strength 
functions for the simplest charge-exchange mode: the isobaric analog resonance 
(IAR) J n = + . The one-body Fermi transition operator reads: 

?£ = I>=f(0- (ii) 

i=l 

For N > Z nuclei TT- simply changes a neutron into a proton without spin- 
flip or change in orbital angular momentum. In Fig. 3 we plot the calculated 
PN-QRPA response to the operator Eq. (11) for 48 Ca, 90 Zr and 208 Pb. The 
strength distributions are dominated by a single IAR peak which corresponds 
to a coherent superposition of proton-particle - neutron-hole excitations. The 
calculated IAR excitation energies (evaluated with respect to the calculated 
ground state of the parent nucleus) are compared with the corresponding 
experimental values from (p,n) charge-exchange scattering data for 48 Ca (7.2 
MeV) [26], 90 Zr (12.0 MeV) [27], and 208 Pb (18.8 MeV)[28]. The agreement 
between the PN-QRPA results and experimental data is very good. We have 
also verified that the calculated strength distributions exhaust the Fermi sum 
rule 

S$- = E\^f\ T p-\^)\ 2 = m\TM) =N-Z, (12) 
/ 

to better than 99.75% (see Fig. 3). 

As a second example we consider open-shell nuclei and plot, in Fig. 4, the 
calculated IAR excitation energies for the sequence of even-even Sn target 
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nuclei with A = 108 — 132. The result of our self-consistent RHB plus PN- 
QRPA calculation are shown in comparison with experimental data obtained 
in a systematic study of the ( 3 He,t) charge-exchange reaction over the entire 
range of stable Sn isotopes [32], and with the relativistic PN-QRPA analysis of 
Ref. [10], in which the phenomenological density-dependent meson-exchange 
interaction DD-ME1 [33] was used both in the RHB calculation of the ground 
states and in the particle- hole (ph) channel of the isovector residual interac- 
tion. Both models reproduce the empirical mass dependence of the IAR, and 
the calculated excitation energies are in very good agreement with available 
data. For DD-ME1 the largest difference between the theoretical and experi- 
mental IAR excitation energies is ~ 200 keV, and an even better agreement 
with data is obtained in the present calculation. This is because the FKVW 
effective interaction reproduces the ground states of Sn isotopes (e.g. binding 
energies) better than DD-ME1 [1]. We note that, because the IAR is a non 
spin-flip mode, the one-pion exchange Eq. (3) and the Landau-Migdal term 
Eq. (4) do not contribute to the matrix elements of the residual interaction 
Eq. (5). The only contribution comes from the isovector channel of the FKVW 
interaction (i.e. the Lagrangian density Eq.(2)). The density- dependent cou- 
pling strengths G^s an d G^v have been adjusted to the (asymmetric) nuclear 
matter equation of state and ground-state properties of finite nuclei. 

We would also like to emphasize a result that was already discussed in Ref. [10], 
namely that sharp, non-fragmented IAR peaks in open-shell nuclei are only 
obtained when the T = 1 pairing interaction is consistently included both in 
the RHB calculation of the nuclear ground-state and in the proton-neutron 
residual interaction. This is simply a consequence of the fact that the pair- 
ing interaction is isospin invariant, and therefore commutes with the Fermi 
operator Eq. (11). 

We note that the Nolen-Schiffer anomaly [29] is absent in the present calcu- 
lation because the attractive Coulomb exchange term is not included in the 
energy density functional that we use in the calculation of ground state proper- 
ties [30]. Further refinements should also consider charge symmetry breaking 
contributions, e.g. from iterated and irreducible pion-photon exchange pro- 
cesses [31]. 

3.2 Gamow-Teller excitations 

The calculated Gamow-Teller (J w = 1 + ) strength distributions for 48 Ca, 90 Zr 
and 208 Pb are shown in Fig. 5. The one-body Gamow-Teller operator reads: 

^ = ES(i)r T (t). (13) 

i=l 
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The corresponding integrated strengths satisfy the Ikeda sum rule: 

S G p I - 1$? = £ \&f\T$- T m\ 2 - T,\&f\T?m\ 2 = 3(iV - Z) . (14) 
/ / 

In addition to the high-energy GT resonance - a collective superposition of 
direct spin- flip (j = I + ~ — > j = I — |) transitions - the response functions 
display a concentration of strength in the low-energy tail. The transitions 
in the low-energy region correspond to core-polarization (j — I ± ~ — > j — 
I ± 2), and back spin-flip (j = I — ~ —> j = I + |) neutron- hole - proton- 
particle excitations. The calculated GTR are compared with the experimental 
excitation energies: 10.5 MeV for 48 Ca [26], 15.6 MeV for 90 Zr [27,28], and 
19.2 MeV for 208 Pb [15,16,17]. Although one of the parameters of the Landau- 
Migdal interaction has been adjusted to reproduce the GTR excitation energy 
in 208 Pb, we find a very good agreement with experiment also for 48 Ca and 
90 Zr. The integrated strengths satisfy the Ikeda sum rule with high accuracy. 
This is an important test of the internal consistency of our relativistic PN- 
RPA. We note that the Ikeda sum rule is exhausted by the calculated GT 
strength only when the relativistic RPA/QRPA space includes both the ph 
excitations formed from ground-state configurations of the fully or partially 
occupied states of positive energy, and the empty negative-energy states from 
the Dirac sea [10]. The contribution of these configurations to the Ikeda sum 
rule is of the order of 8 — 10%. 

Finally, for the sequence of even-even Sn target nuclei, we compare in Fig. 6 
the PN-QRPA predictions for the GTR excitation energies with experimental 
data from Sn( 3 He,t)Sb charge-exchange reactions [32], and with the results 
obtained with the DD-ME1 meson-exchange effective interaction in Ref. [10]. 
The same T = 1 Eq. (6) and T = Eq. (7) pairing interactions have been 
used in both PN-QRPA calculations. Both models reproduce the isotopic trend 
of GTR excitation energies. For the individual nuclei the level of agreement 
with data varies. Below A = 120 the GTR calculated with DD-ME1 plus 
the residual interaction of Ref. [10] are closer to the experimental energies, 
whereas for A > 120 the present PN-QRPA calculation using the FKVW 
parameterization plus the residual interaction Eq. (5) predicts GTR in better 
agreement with data. 



4 Conclusions 

In summary, we have extended our approach to nuclear density functional 
theory, based on chiral dynamics and the symmetry breaking pattern of low- 
energy QCD, to the description of charge-exchange excitations in finite nuclei. 
Starting from the relativistic nuclear density functional (FKVW) introduced 
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in Ref. [1], the proton-neutron (quasiparticle) random phase approximation 
has been formulated and applied in a study of the role of chiral pion-nucleon 
dynamics in excitation modes which include spin and isospin degrees of free- 
dom. In addition to the isovector channel of the FKVW effective interaction, 
the RPA residual interaction includes the direct one-pion exchange and the 
zero-range Landau-Migdal term. 

The density dependence of the Landau-Migdal parameter g' (kf) which charac- 
terizes the spin-isospin channel of the isotropic part of the quasi-nucleon inter- 
action at the Fermi surface, has been calculated in Ref. [12] using in-medium 
ChPT. In the present work one of the short-distance constants has been fine- 
tuned in such a way that the relativistic PN-QRPA calculation reproduces 
the excitation energy of the Gamow- Teller resonance (GTR) in 208 Pb. Both 
approaches lead to very similar values of g' (kf). In addition to the residual 
interaction in the ph channel, the model includes both the T = 1 and T = 
pairing channels. PN-QRPA calculations have been performed for the J 71 = + 
and J n = 1 + charge exchange modes in 48 Ca, 90 Zr and 208 Pb, and for a se- 
quence of even-even 108 - 132 g n nuclei. Results for the excitation energies of IAR 
and GTR, and especially for the isotopic trend in the chain Sn nuclei, are in 
very good agreement with available data and with calculations performed us- 
ing modern phenomenological relativistic density-dependent meson-exchange 
effective interactions. 

This analysis has shown that a nuclear energy density functional based on 
chiral effective field theory provides a consistent microscopic framework not 
only for ground-state properties, but also for complex excitations of the nuclear 
many-body system. The specific density dependence of spin-isospin dependent 
particle-hole interactions, as derived from in-medium chiral dynamics, gives a 
successful description of both isobaric analogue and Gamow- Teller modes. 
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Fig. 1. Left panel: the Landau-Migdal parameter g' as function of the Fermi mo- 
mentum kf for several values of b§ T in the interval [—5.1,-9.5]. In the panel on the 
right the corresponding Gamow- Teller strength functions in 208 Pb are shown. The 
experimental position of the GTR is reproduced with the choice: b^ T = —7.30. 
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Fig. 2. The dimensionfull Landau-Migdal parameter g' as function of the Fermi 
momentum kf (solid curve). The values around saturation point are compared 
with those reported in Refs. [18] (dark gray rectangle), [19] (black rectangle), 
and [20] (light gray rectangle). The black dots are the values obtained with the 
Brueckner-Hartree-Fock calculation of Ref. [21]. The dimensionless G' values from 
Refs. [18,19,20,21] have been divided by the density of states at the Fermi surface 
N = 2TT- 2 k f0 M* rj 1 fm" 2 . 
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Fig. 3. PN-RRPA J 71 " = + strength distributions. The excitations of the isobaric 
analog resonances are compared with experimental data (arrows) for 48 Ca [26], 
90 Zr [27], and 208 Pb [15,28]. The integrated strengths are compared to the mod- 
el-independent Fermi sum rule Eq. (12). 
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Fig. 4. RHB plus PN-QRPA results for the isobaric analog resonances of the se- 
quence of even-even Sn target nuclei, calculated with the FKVW + Gogny effective 
interaction and compared with experimental data [32], and with the calculation of 
Ref. [10]. 
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Fig. 5. Gamow-Teller strength distributions for 48 Ca, 90 Zr and 208 Pb. PN-RQPA 
results are shown in comparison with experimental data (arrows) for the GTR 
excitation energies in 48 Ca [26], 90 Zr [27,28], and 208 Pb [15,16,17]. 
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Fig. 6. Excitation energies of Gamow- Teller resonances for the sequence of even-even 
108-I32g n target nuclei. The present relativistic PN-QRPA calculation is compared 
with the results of Ref. [10], and with experimental data from Ref. [32]. 
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